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Abstract 

The paper addresses parametric inequality systems described by polynomial functions in finite dimen¬ 
sions, where state-dependent infinite parameter sets are given by finitely many polynomial inequalities 
and equalities. Such systems can be viewed, in particular, as solution sets to problems of generalized 
semi-infinite programming with polynomial data. Exploiting the imposed polynomial structure together 
with powerful tools of variational analysis and semialgebraic geometry, we establish a far-going extension 
of the Lojasiewicz gradient inequality to the general nonsmooth class of supremum marginal functions as 
well as higher-order (Holder type) local error bounds results with explicitly calculated exponents. The 
obtained results are applied to higher-order quantitative stability analysis for various classes of opti¬ 
mization problems including generalized semi-infinite programming with polynomial data, optimization 
of real polynomials under polynomial matrix inequality constraints, and polynomial second-order cone 
programming. Other applications provide explicit convergence rate estimates for the cyclic projection 
algorithm to find common points of convex sets described by matrix polynomial inequalities and for the 
asymptotic convergence of trajectories of subgradient dynamical systems in semialgebraic settings. 
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1 Introduction 

This paper is largely devoted to polynomial semi-infinite optimization and related topics being revolved 
around deriving explicit error bounds for infinite parametric inequality systems with (real) polynomial data 
as well as their various applications to stability analysis in optimization and convergence of algorithms. 
The imposed polynomial structure allows us to widely use powerful tools of semialgebraic geometry, while 
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parametric inequalities naturally call for applying constructions and results of variational analysis and gen¬ 
eralized differentiation. Needless to say that the seminal contributions by Terry Rockafellar to variational 
analysis and optimization are difficult to overstate, and it is our honor to dedicate this paper to him. 

The primary attention of this paper is paid to the parametric inequality systems 

S := {x G MJ^\ fi{x,y) < 0 for all y G Y{x), I = 1,L} , (1.1) 


where each function fi: K” x R.™ as I G L for the given natural number L G N is a polynomial, and where 
Y : R." R™ is a set-valued mapping that is also described by finitely many polynomials via inequality and 
equality constraints. Systems of type (HH) naturally arise as feasible solution sets in problems of generalized 
semi-infinite programming, second-order cone programming, robust optimization, and matrix inequalities 
with polynomial data; see below for more details and applications. 

One of the most important issues associated with the inequality systems (HID is establishing the so-called 
error bounds. Given x G R", recall that a (local) error bound of S with a Holder exponent r G (0,1] holds 
for f at x if there exist constants c > 0 and e > 0 such that 


dist(a:, S) < c 


sup 

yev (a:). 
1<Z<L 


fiix,y) 


J + 


for all X with lO — x\\ < s, 


( 1 . 2 ) 


where dist(a:,5) signifies the Euclidean distance between x and 5, and where [a]+ := niax{a,0}. The 
supremum in dm is obviously achieved and it can be replaced by ‘max’ if Y(x) is closed and bounded. 

The study of error bounds has attracted a lot of attention of many researchers over the years and has 
found numerous applications to, in particular, sensitivity analysis for various problems of mathematical 
programming, termination criteria for numerical algorithms, etc. We refer the reader to [32] for an excellent 
survey in these directions and to the more recent papers [H HU dll El EH EH [36] with the bibliographies 
therein. It is worth noting that the major attention in the aforementioned and many other publications on 
error bounds has drawn to the case of linear rate {j = 1), where this issue is related to metric regularity and 
subregularity notions in basic variational analysis. Our main interest in this paper concerns fractional/root 
error bounds 0 < r < 1 in dm- For the case of finite and fixed sets Y{x) = Q in (II. 2p . some results in this 
direction have been obtained in, e.g., [g dni EQl EHESl Ea ED] with various applications therein. 

It is proved [25] in this finite case of Y(x) = O, by using the cerebrated Lojasiewicz gradient inequality 
[24] . that (|1.2p holds with some unknown exponent t G (0,1) for polynomial systems (|l.lll . Employing 
advanced techniques of variational analysis, we have recently derived in this case [22] several error bounds 
with exponents explicitly determined by the dimension of the underlying space and the number/degree of 
the involved polynomials. The techniques and results developed in [22] allowed us to resolve several open 
questions raised in the literature, which include establishing explicit Holder error bounds for nonconvex 
quadratic systems and higher-order semismoothness of the maximum eigenvalue for symmetric tensors. 

The primary goal of this paper is to obtain explicit error bounds of type m for polynomial inequality 
systems with infinite and variable sets Y{x) C R"*. Besides undoubted importance of these issues for their 
own sake, we have been motivated by applications of infinite polynomial systems and error bounds for them 
to higher-order stability and convergence rates of algorithms in optimization-related areas as well as in 
asymptotic analysis of dynamical systems, where estimates of type (HID with infinite sets Y{x) are crucial. 

As the reader can see below, deriving error bounds for the case of infinite and variable sets Y {x) in (11.11) is 
signihcantly more involved in comparison with our developments for finite systems in [22] . First we present 
the following three-dimensional example showing that the error bound p.2ll may fail for any r > 0 for infinite 
polynomial inequality systems even in the case of constant sets Y{x) = H. 

Example 1.1 (failure of Holder error bounds for infinite polynomial systems). Consider the 
polynomial system of type HID containing only one inequality given by / : R x R^ —R in the form 


f{x,y)<0 with y G H, where f{x,y):=xyi+y 2 ,y = {yi,y 2 ), 
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where the infinite set fl is constructed as follows. Take the C°“-smooth function of one variable 


(l){x) 


e ^ if X ^ 0, 

0 if a; = 0 


and define the set C by the conditions 


^ ■= {{yi,y2)\^t e [-0.5,0.5] 


such that 2/1 = ip'it) and 2/2 




(1.3) 


where </>' stands for the derivative of </>. Since </> is C°°-sniooth, the set 17 in (11.31) is nonempty and compact 
as the image of a compact interval under a continuous mapping. We claim that 


ma.xf{x,y) = (j){x) whenever a: G [—0.5, 0.5]. (1.4) 

yen 

Indeed, for any x £ [—0.5,0.5] and y G 17 there is t G [—0.5, 0.5] with yi = 4>'{t) and 1/2 = Thus 

f{x, y) = x(j)'{t) + - t(j)'{t) = (j){t) + (a; - t)<(i'(t). 


Applying the second-order Taylor expansion to the function cj) tells us that 

(/>(a;) = (l){t) + (j) {t){x — t) (x — ty, for some a G [a;, 7] C [—0.5, 0.5]. 

Note that = (4a“® — 6a“"‘)e“^ = (2 — 3a^) > 0, where the last inequality holds due to the 

choice of a G [—0.5,0.5]. Hence we get the relationships 

/(x, y) = cj)(t) -|- (x — t)4>'(t) < (j){x) for all x G [—0.5,0.5], y G fl, 


which imply the inequality maxygn/(x, y) < <(>(x) whenever x G [—0.5,0.5]. On the other hand, it follows 
from the above constructions of / and 17 that 

max/(x, y) > x(()'(x) -b ((()(x) — X(()'(x)) = (()(x) for each x G [—0.5, 0.5], 
yen 

which therefore justifies the claim in (11.41) . Having this in mind, consider the set 

S' := {x G M| /(x, y) < 0 for all y G 17} (1.5) 


and observe that S D [—0.5,0.5] = {xj cj){x) < 0} = {0}. Let us now check that the local Holder error bound 


(11.21) fails for (11.51) at x = 0 for any exponent r > 0. To see it, take Xk '■= 


ydn(fe+T) 

for all large k that d (xfc, S) = d (xfc, S D [—0.5,0.5]) = Jx^J. This allows us to conclude that 


0 as fc G N and get 


d{xk,S) 

[ma.Xy(zQf{xk,y)]l 


\xk\ 


{k + l)T 

— , = —>■ 00 as fc —>■ 00 

VHk + 1 ) 


whenever r > 0 is chosen, and thus the error bound fails for the system S in dEni). 

In what follows we prove that such a situation does not emerge if the sets Y (x) in (11.11) are described by 


Y{x) = {y G g^{x,y) < 0, i = 1,... ,r, and hj{x,y) =0, j = 1,... ,s}, (1.6) 

where gi and hj are polynomials. It is shown in Section 4 that r in (|1.2I) is explicitly calculated in terms 
of degrees of the polynomials and dimensions of the spaces in question. The key of our analysis is a new 
nonsmooth extension of the Lojasiewicz inequality to the class of supremum marginal functions 


(j){x) := sup /(x,y) 

yeY(x) 
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described by polynomials / and gi , hj in (HH). This is done in Section 3 by using powerful tools of variational 
analysis and semialgebraic geometry reviewed in Section 2. 

Sections 5 and 6 are devoted to applications. In Section 5 we develop quantitative higher-order stability 
analysis for remarkable classes of polynomial optimization problems: generalized semi-infinite programming, 
optimization of matrix inequalities, and second-order cone programming. Finally, Section 6 contains explicit 
estimates of convergence rates for the cyclic projection algorithm to solve feasibility problems for convex 
sets described by matrix polynomial inequalities and also for asymptotic analysis of subgradient dynamical 
systems governed by maximum functions with polynomial data. 


2 Tools of Variational Analysis and Semialgebraic Geometry 


This section briefly discusses some tools of generalized differentiation in variational analysis and of semial¬ 
gebraic geometry widely used in the paper. Throughout this work we deal with finite-dimensional Euclidean 
spaces labeled as R" and endowed by the inner product {x, y) = x'^y. The symbol ]Be(a:) (resp. Be(a;)) stands 
for the open (resp. closed) ball with center x and radius e > 0 while B (resp. B) stands for the open (resp 
closed) unit ball centered at the origin in R". Given a set A C R”, its interior (resp. boundary, convex hull, 
and conic convex hull) is denoted by int A (resp. bdA, coA, and cone A). 

Starting with variational analysis, recall first two subdifferential notions needed in what follows. The 
reader can find more information and references in the books EZllsi]. 

Given a function /: R" —> R continuous around x, the proximal subdifferential of / at a: is 


dpf{x) u G 


^ f{x + h)- f{x)- {v,h) 

hminf -—- > —oo 

||?i|Ko. MO \\h\\^ 


( 2 . 1 ) 


The limiting subdifferential of / at a; (known also as the general, basic or Mordukhovich subdifferential) is 


df{x) := {u G R"| —>■ a:, Vk ^ v with Vk G dpf{xk), k G N}. (2.2) 

We clearly have dpf{x) C df{x), where the first set may often be empty (while not so in a dense sense), 
but the second one is nonempty for any locally Lipschitzian function. Furthermore, the set dpf{x) is always 
convex but may not be closed, while df{x) is closed but may often be nonconvex. Both subdifferentials (12.11) 
and (12.2|) reduce to the gradient {V/(a:)} for smooth functions and to the subdifferential of convex analysis 
for convex ones. A significant advantage of the limiting subdifferential (12.2|) is full calculus in the general 
nonconvex setting that is based on variational and extremal principles; see [HIM] for more details. 

The major variational notion used below is the limiting subdifferential slope oi f at x defined via (12.21) by 

mf{x) :=mf {\\v\\\v € df{x)}, (2.3) 

where inf{0} := oo. It reduces to the classical gradient slope Tn/(a:) = ||V/(a:)|| for smooth functions. 

Let us next formulate some continuity notions for set-valued mappings E : R" M R”^; see, e.g., [53]. It 
is said that F is outer semicontinuous (o.s.c.) at x if for any sequence {xk,yk) converging to (x,ff) with 
yk G F{xk) we have y G F{x). Further, F is inner semicontinuous (i.s.c.) at x if for any sequence Xk ^ x 
and any y G F{x) there are yk G F{xk) as fc G N satisfying yk —>■ y. We say that F is o.s.c. or i.s.c. around 
X if it has this property at every a; in a neighborhood of x. 

Finally, we present some notions and facts from (real) semialgebraic geometry following [5]. It is said that: 
• Ac R" is a semialgebraic set if it is a finite union of subsets given by 


{a; G R"| fi{x) =0, i = l,...,q, and ffx) >0, i = g-k 1,... ,p}. 
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where all the functions ft, i = 1,... ,p, are polynomials of some degrees. 

• /: A —)> i? is a semialgebraic mapping if it maps one semialgebraic set A C K" to another one B C K™ 
and its graph {(x, y) ^ Ax B\ y = /(x)} is a semialgebraic subset of M" x R™. We say that f: A ^ B 
is locally semialgebraic around x if there exists a neighborhood V of the point (x,/(x)) such that the set 
V n {(x, y) G A X B\ y = /(x)} is semialgebraic in R" x R™. 

The class of semialgebraic sets is closed under taking finite intersections, finite unions, and complements; 
furthermore, a Cartesian product of semialgebraic sets is semialgebraic. A major fact concerning the class 
of semialgebraic sets is given by the following seminal result of semialgebraic geometry. 

Tarski-Seidenberg Theorem. Images of semialgebraic sets under semialgebraic maps are semialgebraic. 

We also need another fundamental result taken from [U Theorem 4.2], which provides an exponent estimate 
in the classical Lojasiewicz gradient inequality for polynomials. For brevity we label it as: 

Lojasiewicz Gradient Inequality. Let f be a polynomial on R" with degree d G N. Suppose that /(O) = 0 
and V/(0) = 0. Then there exist constants c,e > 0 such that for all x G R" with ||x|| < e we have 

||V/(x)|| > c|/(x)|i-^, where T = Agin, d)-^ and .^(n, d) ^ d > 2^ 

3 Lojasiewicz Inequality for Supremum Marginal Functions 

The main aim of this section is extending the Lojasiewicz gradient inequality (12.411 to the following class of 
(polynomial) supremum marginal functions given by 

:= sup /(x,y), (3.1) 

VGY(x) 

where /: R" x R"* —>■ R is a polynomial of degree at most d, and where the set-valued mapping Y: R" R™ 
is defined in (|1.6|) by polynomials gi: R" x R™ —R for f = 1,..., r and hj : R” x R™ —>■ R for j = 1,..., s 
of degree at most d. Functions of type dSU) are intrinsically nonsmooth, and thus deriving a nonsmooth 
counterpart of (12.41) for them requires the usage of an appropriate subdifferential of (j). The reader will see 
below that a nonsmooth version of the Lojasiewicz inequality for (f in terms of the limiting subdifferential 
slope m 0 (x) from ( 12 .3L which replaces the gradient norm in ( 12 .4|) . plays a key role in establishing our 
Holder-type local error bounds and their subsequent applications in this paper. 

To proceed in this direction, we have to calculate the limiting subdifferential of (j) from (13.11) in terms of 
its initial data, which is not an easy task by taking into account that the sets Y (x) are infinite and variable. 
When T(x) = LI, (|3.1I) reduces to the supremum function for which the most recent subdifferential results 
can be found in [28] ; see also the references therein. When the supremum in ()3.ip is replaced by the infimum, 
we arrive at the class of (infimum) marginal functions well investigated in variational analysis [27] . Needless 
to say that the supremum and infimum operations are essentially different in unilateral analysis and that 
(lower) subdifferential properties under supremum are significantly more involved and challenging. 

We first show that the functions of type (13.11) enjoy the following important properties. 

Proposition 3.1 (supremum marginal functions are locally semialgebraic and upper semicon- 
tinuous). Let the sets T(x) in ()1.6I) be nonempty and uniformly bounded around some point x G R”. Then 
the function </> in dSI]) is locally semialgebraic and locally upper semicontinuous (u.s.c.) around x. 

Proof It follows from the Tarski-Seidenberg Theorem and the assumptions made that 4> is well defined and 
semialgebraic around x. Let us verify that cj) is u.s.c. around x. Select a number e > 0 and a compact set 
C C R”^ so that 0 ^ F(x) C C for all x G ]B£(x). We claim that if is locally u.s.c. on Be(x). Indeed, assume on 
the contrary that </> is not u.s.c. at some x G ®£(x) and then get a sequence x^ —>■ x and a number d > 0 such 
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that (j){xk) > (j)[x)-\-5. Using the continuity of /, the closedness of Y{x) for all x, and the inclusion Y(xt) C C 
allows us to find yk € Y{xk) with 4){xk) = f(xk,yk)- By passing to a subsequence we may suppose without 
loss of generality that yk —t y. It is easy to check that y £ Y{x) due to the continuity oi gi as i = 1,... ,r 
and hj as j = 1,..., s in (HH). Since / is continuous, we have 4>ixk) = f{xk,yk) —t fix,y) < <?!>(a;), which 
contradicts the assumption that 4‘{xk) > 4’i.'^) + Thus 4> is locally u.s.c. around x. □ 

Next we study the lower semicontinuity (l.s.c.) and Holder continuity of the function (j) (13.11) under rather 
mild assumptions. Define the argmaximum set important for our further analysis by 

M{x) := {y £Y{x)\ ^{x) = f{x,y)} (3.2) 

and recall that dist(H; B) := inf{||a — b\\ \ a £ A, b £ B} for A, B C K™. 

Proposition 3.2 (continuity of supremum marginal functions). Given x £ R", suppose that the sets 
Y(x) are nonempty and uniformly bounded around x. The following assertions hold: 

(i) If dist{M{x);Y(x)) —^0 as x ^ x, then (j) is l.s.c. {and hence continuous) at x. Thus the i.s.c. property 
of Y at X ensures that cj) is continuous at this point. 

(ii) If there are numbers e, L > 0 and a £ (0,1] such that 

dist (M(a:i); y(a: 2 )) < T||a;i — X 2 ||“ for all xi,X 2 £V>s{x), (3.3) 

then (j) is Holder continuous on Be(a;) with order a. 

Proof. Since the sets Y(x) are nonempty and uniformly bounded around x, there exist a compact set C and 
a number e > 0 such that Y(x) C C for all x £ ]Be(x). Let us first verify (i) while arguing by contradiction. 
Suppose that dist(M(x);U(x)) —>■ 0 as x —>■ x and that (j) is not l.s.c. at x. Then there exist a constant 
v > 0 and a sequence Xfc —>■ x such that 4>{xk) < 4i{x) — v. Since d{M{x)]Y{xk)) —t 0 as fc —>■ oo, we 
find yk £ Y(xk) and yk £ M(x) with ||j/fe — j/fe|| —t 0. Using the boundedness of {yk}keN passing 
to a subsequence if necessary allow us to find y G R™ with yk ^ y and thus yk —t y. It is easy to see 
that y £ Y(x) and that 4){x) = f{x,yk) —t f{x,y), which yield in turn that y £ M(x). Hence we have 
4>{xk) > f{xk,yk) —t fix^y) = 4>ix) as fc —>■ oo while contradicting the assumption on (f{xk) < 4>{x) — v. 
This verifies that (j) is l.s.c. at x and hence its continuity there by Proposition 13.11 The second statement in 
(i) follows from the first one and the definitions. 

To justify now assertion (ii), suppose that (13.31) holds and pick any xi,X 2 G ]B£(x). Since the sets M{x) 
and Y{x) are compact for all x G ]Be(x), we deduce from (13.311 that there are yi £ M{xi) and j /2 G Y{x 2 ) 
satisfying \\yi — 2 / 2 II < L||a;i — X 2 ||“. Taking into account the Lipschitz continuity of the polynomial / on 
B£(x) X C with some constant f > 0, it follows that 

(t){xi) - (^(x 2 ) < fixi,yi) - /(x 2 ,j/ 2 ) < ii\\xi - X 2 II + ||yi - 2 / 2 II) 

< f(||xi - a; 2 || + L||xi - X 2 |r) = f((2e)^““ + L)||xi - X 2 |r. 

Similarly we conclude that 4>{x2) — (f{xi) < f((2e)^““ + T)||xi — X 2 ||“, which verifies the Holder continuity 
of (j) on B£(x) with order a and thus completes the proof of the proposition. □ 

Remark 3.3 (effective conditions for continuity of marginal functions). Proposition I3.2f il tells us 
that (f in dSlD is continuous at x provided that Y is i.s.c. at this point, which surely holds if Y is locally 
Lipschitzian around x in the standard/Hausdorff sense with taking into account that the sets Y{x) from (11.61) 
are closed, nonempty, and uniformly bounded. Necessary and sufficient conditions for the local Lipschitz 
property of Y around x are given in m Corollary 4.39] in terms of the gradients of the constraint functions 
gi,hj at (x,y) for any y £ Y{x). In this way we get effective conditions ensuring the continuity of (j) at (in 
fact around) x via the initial data of (11.61) . On the other hand, |271 Theorem 3.38(iv)] justifies the Lipschitz 
continuity of ft around x under the i.s.c. property of the argmaximum mapping (13.21) around (x, y) with some 
fixed y £ M{x) and the Lipschitz-like (Aubin, pseudo-Lipschitz) property of Y around this pair (x, y), which 
is effectively characterized in HZl Corollary 4.39] via the initial data. 
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To proceed further, we need the following qualification condition imposed at a reference point relative to 
some (in fact optimal) subset of the constraint set Y(x) in (13.ip . 

Definition 3.4 (marginal constrained qualification). Given x G K" and Y(x) from (11.61) . we say that 
the MARGINAL Mangasarian-Fromovitz CONSTRAINT QUALIFICATION (MMFCQ) holds at X relative to 
some subset fl <zY(x) if there is a vector f G R" such that 

r 1 

r s 

'^XiWxgz{x,y)+ '^KjVxhj{x,y) £, > 0 for any y G (3.4) 

_i=i j=i 

whenever Lagrange multipliers (A, k) G ]R!j_ x R® \ {0} satisfy the conditions 


'^XiVyg^{x,y)+ '^KjVyhj{x,y) = 0 and X^g^{x,y) = 0, i = l,...,r. (3.5) 

i=i j=i 

Remark 3.5 (relationships of MMFCQ with EMFCQ and MFCQ). The defined MMFCQ for the 
marginal supremum functions is motivated by the extended Mangasarian-Fromovitz constraint qualifica¬ 
tion (EMFCQ) introduced in [T7] for generalized semi-infinite programs (GSIPs) discussed below in Section 5. 
Reformulating EMFCQ for (13.11) . we see that it involves the cost function / and also requires the validity 
of dSH) for a larger set of Lagrange multipliers in comparison with MMFCQ; so MMFCQ is weaker than 
EMFCQ in general. Observe also that when the functions gi and hj do not depend on x, the formulated 
MMFCQ condition means that there is no (A, k) G R!j_ x R® \ {0} satisfying (13.51) . which is equivalent to the 
conventional Mangasarian-Fromovitz constraint qualification (MFCQ) on Q in the sense that 

r s 

XiVg^iy) + Y = 0 : 

i=l j=l 

Xt > 0, Xigiiy) =0, z = 1,..., r, and /Cj G R, j = 1,..., s 

for any y G LI. As we see in the subsequent sections, the imposed MMFCQ holds automatically for important 
classes of polynomial optimization and related problems arising in applications. 

It is worth emphasizing that the majority of our applications below requires the validity of MMFCQ for 
the case when LI in Definition 13.41 is chosen as the argmaximum set (13.21) at the reference point. Consider 
further the standard Lagrangian function for dSH) with the negative sign of 7 taken due to the maximization 


Xi = 0 , Kj = 0 


(3.6) 


r s 

■X^’{x,y,'y,X,K) := -7/(1, y) + ^ A,g,(a;, y) + ^ «;j/ij(a;, y). 

i=i j=i 


(3.7) 


Observe that any y G M (x) is a minimizer of the following nonlinear program: 


min—/(x, •) subject to yi(x, •) < 0, z = 1 ,..., r, and hj{x, ■) = 0, j = 1,..., s. (3-8) 


Applying the classical Lagrange multiplier rule in the Fritz John form to (|3.8I) tells us that the set FJ{x,y) 
of multipliers ( 7 , A, k) G R+ x R!j_ x R® satisfying the conditions 

r s 

7 + + X! ^ Aigi(x,y) = 0, i = 1,... ,r, (3.9) 

i=i j=i 

is always nonempty. Given a G [0, 00 ], consider also the set FJa{x, y) of ( 7 , A, n) G R+ x R!,. x R® satisfying 


7 + ^ Ai + ^ |kj| < a, Vyif(x,y, 7 , A, At) = 0, and Aiyi(x, y) = 0, i = 1,..., r. (3.10) 

i=l j=l 

Using the Lagrangian description, we now show that the MMFCQ property relative to O = M{x) from 
(13.21) is robust in the general setting of Proposition 13.2r iL 
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Proposition 3.6 (robustness of MMFCQ). Given x S R", assume that the sets Y(x) are nonempty and 
uniformly bounded around x. If MMFCQ holds at x relative to M[x) and if 4> is l.s.c. at x, then there exists 
<5 > 0 such that MMFCQ is satisfied at any point x G ® 5 ( 2 :) relative to M(x). 

Proof. Proposition 13.11 tells us that (j) is continuous at x. Supposing on the contrary that MMFCQ is 
not robust, i.e., there are —>■ a; as fc —>■ oo such that MMFCQ fails at Xk relative to M{xk). Take 
any ^ G R" and find yu G M{xk) and {Xk,Hk) € R!j_ x R® \ {0} satisfying '\/yF!f{xk,yk,0,Xk,Kk) = 0, 
{Xk)igi{xk,yk) = 0 for i = 1 ,... ,r, and Vx-^{xk,yk-, 0, Xk,Kk)'^f, < 0. Normalization gives us ||(Afe, Kk)\\ = 1 
for all fc G N. Using now the uniform boundedness of F(a:) around x, we select subsequences yk ^ y G Y{x) 
and (Afe, Kk) —>■ (A, R) with ||(A, /t)|| = 1. It follows from the continuity oi (f at x and from 4){xk) = f(xk,yk) 
that (fix) = f{x,y), i.e., y G Mix). Furthermore, passing to limit as fc —>■ oo yields Xigiix,y) = 0 for 
i = 1,..., r, Vy.if (a;, y, 0, A, R) = 0, and (x, y, 0, A, R)^£, < 0. Since ^ G R" was chosen arbitrarily, this 
contradicts the assumed MMFCQ at x relative to M{x) and thus completes the proof. □ 

The next result important for its own sake plays a crucial role in deriving the extended Lojasiewicz 
inequality for supremum marginal functions. It explicitly evaluates the limiting subdifferential of for such 
functions via the initial data and the corresponding Lagrange multiplies being significantly different from 
the preceding results and techniques of |28] even for the constant mapping Y in (13.11) considered therein. 
Theorem 3.7 (limiting subgradients of supremum marginal functions). Given x G R", suppose 
that the sets U(x) in (11.61) are nonempty and uniformly bounded around x, that MMFCQ holds at x relative 
to Mix), and that cf is l.s.c. around x. The following assertions hold: 

(i) There exists £ > 0 such that for any x G BE(a;) and v G dcfix) we can find G Mix) and 

G FJaoix,y^'^^), q = I,...,n + I, satisfying the conditions 

n+1 n+1 

= and ^ 7 ( 9 )=!, (3.II) 

q=l q=l 

where and F Joa are defined in (13.711 and (13.101) . respectively. 

(ii) Given L > 0, there are positive numbers s,a such that for any x G ]Be(x) and v G dcfix) D (LB) we can 
find y^^^ G Mix) and ( 7 ^'^^ A^®\ G LJq(x, q = I,...,n + I, satisfying (13.III) . 

Proof. It follows from Proposition 13.61 that there is £ > 0 such that MMFCQ holds at any point x G Beix) 
relative to Mix). Moreover, Proposition 13. 1 1 and the l.s.c. assumption on (f ensure that this function is con¬ 
tinuous on Be(x). Let us first evaluate the proximal subdifferential (12.11) and then the limiting subdifferential 
(12.21) of the supremum marginal function (13.11) at each fixed x G B^)!:). 

Claim 1: For any proximal subgradient v G dpcfix) we have 

(u,-1) G cone{(Va;.if(x,y, 7 , A, k),- 7 )| y G M(x), ( 7 , A, k) G LJ(x,y)} C R"'+^ (3.12) 

To verify (13.1211 . deduce from (12.111 that there are constants 7 > 0 and (5 G (0, £ — ||x — x||) satisfying 

(fiu) — (fix) > iv,u — x) — r]\\u — xW"^ whenever u G B 5 (x). (3.13) 

The assumptions imposed on the mapping Y ensure that (fiu) = max{/(M, y)| y G U(m)} for all u G R", 
which tells us by (13.131) that the pair (x, (fix)) is a local minimizer of the following GSIP: 

min [z — (ri, M — x)-I-ylfy — x|p] subject to fiu,y) — z<0 for all yGYiu). (3-14) 

{u,z) 

Applying to (13.1411 the necessary optimality condition from [T3 Theorem 1.1], we find y G N, y^'^^ G M(x), 
G LJ(x, y^'^^), a > 0, and > 0 for g = 1,... ,y so that 

p p 

a + = 1 and a(—r;, 1) -I- (Va,.5f (x, y^‘^\ 7 *-'^^ A^®\ — 7 ^'^^) = (0, 0). (3.15) 

q=l q=l 
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To justify (|3.12D . it remains to show that a > 0 in f|3.15p . Indeed, assuming the contrary tells us that 


p p 

y]y«) = land = (0,0) (3.16) 

9=1 9=1 

with G FJ{x, for q = 1,... ,p. Consider the set / := {g G {1,... ,p}| O}, which 

is nonempty by J2q=i = 1- follows from (13.161) that = 0 for g = 1,... ,p. Hence = 0 and 

(0, kI'5^1) G FJ{x, for all q& I. Combining this with (13.161) yields 

5 ]/i(«)V,JS^(a:,y(«),0,A(«),s:(«)) = 0, 
q&I 

which contradicts the assumed MMFCQ at x relative to M{x) and thus verifies (j3.12|) for v G dp(j){x). 
Claim 2: Any limiting subgradient v G d4>{x) satisfies (|3.12|1 and thus (|3.1ip . 

Fix V G d4){x). Since f is continuous on ]Be(a:), we find by definition (12.211 sequences Xk ^ x and Vk G dp4){xk) 
with Vk ^ V. Using (I3.12p for Vk and applying the Caratheodory theorem to the conic convex hull in (13.121) 
give us € M{xk), €R+, and € FJ{xk,yl?^) ior q = 1,... ,n + 1 such that 

n+1 n+1 

'>^k = ^lJ-k'^y^-^ixk,yk\lk^,\k\'^k^^ and = 1. (3.17) 

9=1 9=1 

Let us show that the sequence of pk ■= iy^k \ • ■ •, is bounded. Arguing by contradiction, suppose that 

^9=1 and define := yL^P [ Ylpti mIT^] ^ for g = 1,..., n + 1. Since the convergent sequence 

of Vk is bounded as Vk ^ v, we deduce from (I3.17|) that 


n+1 

{Xk , yp'^ , ) ^ 0 and 

9=1 


n+1 

9=1 


0 as fc —>• 00 . 


(3.18) 


It follows from the boundedness of yp '^, , and ) for g = 1,..., n +1 that some subsequences 

of them converge to y^‘^\ and (y^'^i, A^'^i, respectively. Letting /c —>• oo in (13.1811 yields 


n+1 

^(cr,, y(«), A(«) ,+)) 
g-1 


n+1 

0 and = 0 

9=1 


(3.19) 


with (yl'Jl, A^'^i,) G FJ{x,y^'^p, G M{x), and X^g^i = 1- Similarly to Claim 1 we arrive at a 

contradiction with the validity of MMFCQ at x relative to M{x), which therefore shows that the sequence 
of pk = {Pk\- ■ ■, ^ i® bounded. 

By passing to subsequences again, we get that fik ^ fJ- & IR+^, yp'^ —>■ y^‘>^ G M{x), and {'yp\ —t 

(y^'^i, A^'^i, G FJ{x, y('^l) as fc —>■ oo for g = 1,..., n + 1. It follows from (13.1711 that 


n+1 n+1 

w = y A^^i, and y = 1, (3.20) 

9=1 9=1 

which gives us (13.1211 for each v G d(j){x) as stated in Claim 2. 

To complete the proof of (i), define y^®^ := := and := and then observe 

that the inclusion (yl'^l, A^'?^ kI"?!) G FJ{x,y^^p yields (y^*^ G FJaoix,y^‘^P and 


n+1 n+1 

i; = yV,Jf’(x,+),y(«),A(«\ «(«)), yyl?) = 1. 

9=1 9=1 


(3.21) 
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To justify finally (ii), pick any v G d(j){x) fl (LB) with x G B£(iE) and find by Claim 2 and the Caratheodory 
theorem such G M(x), G M+ and G FJ{x, for y = 1,..., n + 1 that (13.201) holds. 

Taking into account that all the elements above depend on e and proceeding as in the proof of Claim 2, we 
deduce that , q = l,...,n + l, are uniformly bounded when e is sufficiently small. This means that there 
are numbers a,e > 0 such that 0 < < a for any satisfying (13.201) whenever x G Be(a:). Defining 

7 ( 9 ) := ^( 9 ) 7 ( 9 ), := y( 9 )A( 9 )^ and := as above and then using ( 7 («), e FJ(x,y(«)), 

we conclude that G FJa(x,y^‘^^) and that v satisfies (13.211) . This clearly ensures the validity 

of (13.111) and thus finishes the proof of the theorem. □ 

Now we are ready to establish a subdifferential extension of the Lojasiewicz Gradient Inequality to supre- 
mum marginal functions with polynomial data in terms of the limiting subdifferential slope (IQ) . 
Theorem 3.8 (Lojasiewicz subgradient inequality for polynomial systems). Given x G K", suppose 
that all the assumptions of Theorem 13.71 are satisfied for the marginal supremum function (j) in (13.11) . Then 
there exist positive constants c and e such that we have the estimate 

> c\(j){x) — (j){x)\^~'^ for all x gM!;{x) with r = .^(2n + (m + r + s)(n + 1), d + 2) (3.22) 

where the constant is taken from (EH). 


Proof. It follows from Proposition 13.11 that is continuous around x. Hence there are z/, ei >0 such that 
\4>{x) — 4>{x)\ < V for all x G B£j(a;). Define L := v^~'^ and F: K" x K™ x R x K’’ x K® —>■ R by 

r s 

F{x,y,-i,p,K) := - 7 /(x,y) + ^y-y*(a;,y) + ^Kj/ij(a;,y) (3.23) 

i=l j=l 


and observe that L is a polynomial of degree not greater than d + 2. Define further another polynomial 
P: R" X R'"(’^+i) X R" X R’'("+i) x R®("+i) —^ R of degree not greater than d + 2 by 


9=1 


) + F{x, y("+i), 1 - X! 7^"^^, 
9=1 




(3.24) 


where a; e R", y = (y^^),..., y("+i)) e R™("+i) with y(«) S R"*, 7 = ( 7 ( 1 ),..., 7 (")) e R” with 7 ( 9 ) S R, 
y = (y(i),...,y("+i)) G R’'(’"+i) with y(«) G R’', and k = ...,«(”+!)) G R®("+i) with G R®. For 

each y = (yi,..., fir) G R"" denote y^^l := (yf,..., y^) G R!j_ and consider the set 


K{s) := {(x, y, 7 , y, k) G R" x R™("+i) x R" x R’'("+i) x R®("+i) 

n 

a; G le(x), 7 G R", ^ 7 ^’^ < 1, ■ ■ ■, 7^”^ > 0, 

y = (y(^\ ... ,y("+^)), y = (y(^\ ... ,y("+^)), k = (k^^^ 

y(«) G M{x), q=l,...,n + l, 

G FJ„(a:,y(‘?^), g = 1,... ,n, 

GFJ«(:r,y(-+!))}, 


,K 


(n+l)'l 


£ G (0, £ 2 ), 


(3.25) 


where positive the numbers a and £2 < £1 depend on L and are taken from Theorem I3.7f iil. It can be 
verified that the set-valued mapping K{-) is o.s.c. at £ = 0. Furthermore, for any (x,y,^, fi^ k) G K{0) we 
get from (13.101) and p.25|) the equalities 


n n 

P{x, y, i,y,n) = - ^ i^'^'^fix, y^"*^) - ((1 - X! 7 ^'^^)/(a;, + fiix) 

q=l g=l 

n n 

= 7^'^V(S) - (1 - X! 7^''^) '('(S) + Hx) = -fi(x) + 4>{x) = 0, 

q=l q=l 
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where the second one follows from € M(x) = M(x) as q = 1,... ,n + 1. Applying the Lojasiewicz 
Gradient Inequality (EH) to the polynomial P in (13.241) gives us positive constants rj and c,, such that 

if \\{x',y',j',y',K') - {x,y,j,n,K)\\ <r] (3.26) 

with T := .^(2n + (m + r + s)(n + 1), d + 2). It is easy to see that K{0) is a compact set, and hence we can 
find p G N, G K{Q), and yj/ > 0 for i/ = I,... ,p so that 


K{Q) c y 

This allows us to find a positive number 5 ensuring the inclusion 

p 

K{0) + ® C y , 7 ^") , . 

Employing the o.s.c. property of K{-) at 5 = 0, we find ^3 G ( 0 ,^ 2 ) so that 

p 

K{e) C X(0) + dB C y , 7 ^"), 

k^l 

for all 5 € [0, ^ 3 )- Denoting c := min{l, , c^^p}, deduce from f|3.26p that 

l|VP(a;,y, 7 ,p,K)|| > c|P(x,y, 7 ,p,/c)|^“^ if (x, y, 7 , p,/c) G Ar(e). 


(3.27) 


Pick now e G (0, £ 3 ) and x G B£(x). If m^{x) > L, then (13.221) holds by the choices of L and c. If m^{x) < L, 
then we take any v G d(j){x) n (LB) ^ 0 and deduce from Theorem 13.7l iii that there are y^'^^ G M{x) and 
«('*)) G FJa{x, y^*^), q = l,...,n + I, satisfying 


n+1 


n+1 


^ V,,Jf(x,y('?), 7 ('?),A(«),K(«)) and ^ 7 ( 9 )=!. 

q=l q=l 


(3.28) 


Define next := i = 1,... ,r, and observe the inclusions 


(7(«),[/i('?)]2,K(‘?^) G FJ«(x,y(‘?^), y = l,...,n + I. 

It follows from (j3.28p and the definition of K{e) in (j3.25|) that {x,y,j,fj,,K) G K{s), and thus (j3.27p holds 
for this quintuple. Differentiating (I3.24p with taking (13.2811 into account implies the equalities 

n+1 

^xP{x,y,^,y,K) = '^\7^^{x,y^‘‘\-f^'^\X^‘^\K’'^'>) = v, 

q=l 

VyP(x,y, 7 ,p,K) = (\7y(i)Ji'{x,y^'^\-f^'^\X^'^\K^'^'^),.. ., Vy(n+i)if(x, y("+^\ 7 ^"+^^ = 0, 

\7^P{x, y, 7 , y, k) = {f{x, y'^^) - f{x, y^"+^^),..., /(x, y^")) - /(x, y^"+^^)) = 0 , 

V^P(x,y, 7 ,y,K) = (2y'f^ gi{x,y^‘^'>)j _ ^ = 0, and V«,P(x, y, 7 , y, k) = (^+(x, y^'J^))^ ^ = 0, 

which ensure by (j3.27ll . (|3.28L and (13.2411 the relationships 

||u|| = ||VP(x,y, 7 ,y,K)|| > c |P(x, y, 7 , y, «;)+"" = c|(/)(x) - ^(x)+^ for any v G d(p(x} n (LB). 


Since m^{x) < L, this clearly verifies (13.221) and completes the proof of the theorem. □ 

Remark 3.9 (sharper estimate for equality polynomial systems). A close look at the proof of 
Theorem 13.81 reveals that the exponent r in (13.221) can sharpen to r = i^(2n + (m + s)(n + 1), d + 1)“^ in 
the case when the polynomial system in (frel) is given by only the equalities D(x) = {y G K"*! hj{x,y) = 
0, j = 1,..., s}. This is due to the fact that the polynomials F in (13.231) and P in (13.241) are of degrees at 
most d + 1 instead of d + 2 as in the general case. 
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4 Error Bounds for Parametric Polynomial Systems 


In this section we derive Holder error bounds with explicit exponents for general polynomial systems of type 
(HU), where the underlying state-dependent sets of parameters Y{x) admit the polynomial description (I1.6I1 . 
We also present important specifications of these estimates in some particular cases. The results obtained are 
far-going extensions of those established in [35] in the case of finite and constant sets Y{x) = H. To proceed, 
we need the following lemma on error bounds for locally Lipschitz functions taken from Corollary 2]. 

Lemma 4.1 (sufficient condition for error bounds of Lipschitzian functions). Let f: R" M. be 

continuous around x G hdSf, where Sf := {x\ f{x) < 0}. Assume that there are numbers c,e > 0 such that 
nr/(a;) > c \ f{x)\^~'^ for all x € R." with ||x — irll < e and x ^ Sf. Then we have 

It S 

dist(a;,5'/) <—[/{2 ;)]_,_ whenever ||a; —x||<-. 

We are now ready to establish the main result of this section. To avoid some technical complexity on 
dealing with the l.s.c. property of supremum marginal functions as discussed in Section 3, we assume in what 
follows that the mapping Y in (11.61) is inner semicontinuous around the point in question, which is always 
the case when Y{x) = ^1 for some fixed set C R™; see, e.g., (14.51) below. 

Theorem 4.2 (Holder error bounds for parametric polynomial systems). Let fi, gi, hj : R" x R™ — 

R with I = 1,..., L, i = 1,... ,r, and j = 1,..., s be polynomial functions of degrees at most d. Suppose that 
the mapping Y from (HU) is i.s.c. with nonempty and uniformly bounded values Y{x) around a given point 
X G R" and that MMFCQ from Definition \‘d.4\ holds at this point relative to the argmaximum set 

M(x) := {y € Y{x)\ fi{x,y) = fi{x,y) for some 1 = 1,. ..,L}. (4.1) 

1<1<L 


Then there are numbers c, e > 0 such that we have the local error bound 


dist(a::, 5") < c ( max fi{x, y) ) for all x G R" with ||x — x|| < e, 

\<1<L 


(4.2) 


where S is taken from (HU; and where the Holder exponent t in (HU is explicitly calculated by 


( I%(2n + {m + r + s){n + l),d + 2) L = 1, 

\ .^(2n-I-(m-I-r-I-s-|-2)(n-I-1), d-|-L-I-1)“^, L>2, 


with the constant defined in (EU- If we suppose in addition that MMFCQ holds at every x G R" relative 
to the whole set of parameters Y{x) and that Y are i.s.c on R", then for any compact set K C R” there is 
a number c > 0 such that we have the following global error bound with the same exponent t: 


dist(x, S) < c 


max fiix, y) 

y€Y (x) 

1<1<L 



whenever x € K. 


(4.4) 


Proof. Consider first the case of L = 1, where (f^x) = sup^gy^^^ fi(x, y) is continuous around x by Proposi¬ 
tions O and l3.2l ib Theorem 13.81 gives us numbers co,eo such that for all x G ]B£o(^) have the estimate 

> co\(j){x) — (j)[x)\^~'^ with T = I%{f2n-\-[m + r-\-s){n + l,d-\-2') 


Let us justify the existence of c, e > 0 for which (HU holds with T from (IQ) as L = l. If i G bdS”, i.e., 
(j){x) = 0, we get by Lemma SU that (14.31) is fulfilled with c = and e = eo/2. Since this inequality is 
automatic for x G int S, it remains to examine the case ol x ^ S. It follows from x ^ S = {x \ 4>{x) < 0} 
that (j){x) > 0 and so {(l){x)y > 0. By the continuity of </> and d(-, S) we find Mi, M 2 , e > 0 such that 


{4>{x)y > Ml and dist(g:. S') < M 2 for all x G Be(x), 
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and hence dist(a;, S) < M 2 < ■^4'{xy for all x from this ball. This verifies (14.21) for L = 1 . 

Next we consider the case of L > 2 and define the Lagrange interpolation polynomials 7 ;: K. —>■ R by 

7iW := n l = l,...,L, 

for which 7 /(/) = 1 as / = 1,..., L and ^i{j) = 0 as j € {1,..., L} with j ^ 1. Define the functions 

L 

f{x,y,t) := '^li{t)fi{x,y), g,{x,y,t) := g^x^y), i = l,...,r, 

1^1 

L 

hj{x,y,t) := hj{x,y), j = 1 ,... , 5 , and hs-yi{x,y,t) := - /) 

and observe that / is a polynomial on R" x R™+^ with degree not greater than d + L — 1. Form now the 
modified constraint and argmaximum sets by, respectively, 

Y{x) := {{y,t) e R™ X R| gyx,y,t) <0, i = l,...,r; hj{x,y,t) = 0, j = 1,..., s + l}, 

M{x) := { {y, t) gY (x) | J{x, y, t) = max f{x, y, t )} 

{y,t)eY{x) 

while noting that Y{x) =Y{x) x {1,..., L} for F(x) from (11.61) and that 

max/i(x,j/)= max f{x,y,t), S = {x G W^\ f{x,y,t) < 0 for all {y, t) G Y (x)}. 
f<T7 {y,t)^Yix) 

It can be directly checked that the mapping Y is i.s.c. around x with the nonempty and uniformly bounded 
values and that the MMFCQ property for the modified system holds at x relative to M{x). Applying 
the result for the case of L = 1 obtained above to the system f,gi,hj with (n,m,r, s,d) replaced by 
(n, TO + 1, r, s + 1, d + L — 1), we arrive at (14.21) with t calculated in (14.3p for L > 2. The last statement of 
the theorem (14.41) follows from (14.21) by employing the standard compactness arguments. □ 

Remark 4.3 (sharper error bounds for equality polynomial systems). Based on the observation in 
Remark 13.91 and following the proof of Theorem 14.21 we see that the the exponent estimate r in (14.21) can be 
sharpen to r = t%{2n + {m + r + s + 2 min{L — 1, l})(n + 1), d + L) if the polynomial system Y (x) in (11.61) 
contains only the equality constraints. 

In the rest of this section we derive specifications of Theorem 14. 2 1 in some particular classes of parametric 
polynomial systems and discuss the possibility of sharpening Holder exponents in certain situations. 

Error bounds for infinite polynomial inequality systems 

Let us examine the case of (EB where the sets Y (x) are independent of x being reduced to the constant set 
H := {y e R”"| giiy) <0, i = 1,... ,r; hj{y) =0, j = 1,.. .,s} (4.5) 

described by polynomials with degrees at most d. In particular, such a case covers important models known 
as polynomial matrix inequalities] see below. Observe also that LI from (14.51) can be treated as an infinite 
index set in problems of semi-infinite programming (SIP); see, e.g., [T3]. However, semialgebraic structures 
of index sets have not been specifically investigated in the usual framework of SIP. 

Following the traditional terminology, we say that Lt is regular if MFCQ for (14.51) holds on LI. Let us 
present examples of two important classes of regular sets described by polynomials. 

• Regularity of discrete sets. Considering the discrete set LI := {I,... ,p}, p S N, we can write it as 

H = {p e R| hiy) = 0} with h{y) := {y - l){y - 2) ■ ■ ■ {y - p). 
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Since we obviously have Vh{y) 0 for all 1 / G this set is regular. 

• Regularity of algebraic manifolds. It follows directly from the definition that the algebraic manifold 


n:= {y G hi{y) = 0,...,hsiy) = O}, s < m, 

is regular provided that rank(V/ii(y),..., V/is(y)) = s for all ?/ G fl. 

As an immediate consequence of Theorem 14.21 we arrive at the following extension of the results in 
from finite to to infinite polynomial inequality systems. 

Corollary 4.4 (Holder error bounds for infinite polynomial systems). Let fi: R" x R™ —^ R, 

1 = 1,..., L, be polynomials of degrees at most d, let LI ^ (h given in (14.51) he bounded and regular, and let S 
be defined in (HID with Y (x) = LI. Then for any x G R" there exist constants c, e > 0 such that 


dist(x, S) < c 


max fi{x,y) 

yen 

l<i<L 



whenever ||x 


x|| < 


(4.6) 


where r = .^(2n + (m + r + s + 2 min{L — 1, l})(n + l),d + L + 1) ^ with the constant taken from (12.41) . 
Furthermore, the error bound (HU) holds globally with some c > 0 on any compact set K C R". 

Remark 4.5 (semialgebracity of index sets is essential for Holder error bounds). As follows from 
Example 11.11 the error bound (14.6p can fail if Lt is not semialgebraic. Indeed, this example shows that the 
Holder error bound does not hold whenever r G (0,1) for / and fl therein. It remains to check that the set 
LI from la is not semialgebraic. Supposing the contrary implies by the Tarski-Seidenberg Theorem that 
the function x 1 —>■ maxygn fix, y) is semialgebraic on [—0.5, 0.5], which is not the case. 

Next we present a consequence of Corollarv l4.4l to polynomial matrix inequalities. Recall that an (m x m) 
polynomial matrix inequality of n variables with degree d is represented in the form P(x) F 0 , where the 
notation M fi. 0 means that — M is a positive semidefinite matrix, and where P : R" —^ 5"* is a mapping 
such that each (i, j)th element of its image (P(x))^^., 1 < i, j < m, is a real polynomial of degree at most d. 
We refer the reader to, e.g., [151116] for various applications of such systems. 

Corollary 4.6 (error bounds for polynomial matrix inequalities without regularity assump¬ 
tions). Fix d G N and consider the solution set 


SpMi ■= {a;| Pix) fi 0} 


described by the above polynomial matrix inequality. Then for a compact set K C R" there is c> Q such that 


distix, SpMl) < C ( Ainax(H(x)) 


whenever x G K, 


where r = .^(2n-|-(m-|-l)(n-|-1), d-|-3) ^ and Amax denotes the maximum eigenvalue of a symmetric matrix. 


Proof Define /(x, y) := Vi (-P(^)) Vj^ which is a real polynomial of degree d -|- 2 on R"+™. We have 


P{x) fi 0 4=^ fix,y) < 0 for all y G LI := {y 


= 1 } 


(4.7) 


Note that the set D 7 ^ 0 in (|4.7I) is bounded and regular. Applying now Corollary 14.41 to the inequality 
system on the right-hand side of (14.71) and replacing d by d -I- 2 with taking into account Remark 13.91 ensure 
that for any compact set K we arrive at the relationships 


dist(x, SpMi) < c 


sup /(x,y) 

VMIyP=i 


Amax(-P(a;)) 


with some constant c > 0. This completes the proof of the corollary. 


□ 
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Discussion on sharpness of exponent estimations 


Finally in this section, we discuss the possibility to sharpen Holder exponents in the established error 
bounds for general infinite polynomial systems in some particular settings. Let us first present an example 
of a polynomial matrix inequality, where the exponent obtained in Corollary 14.61 is not the sharpest one. 

Example 4.7 (sharper error bound for matrix inequalities). Let d be an even number. For any 

X = (xi,..., Xn) define Ai{x) := xf and then, for any i = 2,..., n, 


A,{x) := 


-1 xf 


-1 0 

0 0 


+ 


0 0 

0 -1 


X^ -Xi^l 

Consider the following (2n — 1) x (2n — 1) polynomial matrix inequality system: 


0 1 
1 0 


/ Ai{x) 0 0 

0 A2(x) 0 


\ 

0 


0 


-< 0 


Pix) := 

Vo 0 ... A^x) J 

and get by direct calculations that Spmi ■= {x S K"| P{x) ^ 0} = {0} and 

'^max ~ {'^max(^2(^)) I l| 

'-(l+Xi-l) /(I - Xi-l)'^ 


max I max | - 


+ 2 < i < n| 


= max < max ■ 


xf‘^ - Xi-i 


(l+Xj-i) 

2 


{l-Xj-iY , 2d 
4 'X 




2 < i < n 


Observe further that for the parametric family x(t) := ,..., <), t G (0,1), we have 

d{x{t), SpMi) = 0{t) and X-max{P{x{t))) = . Thus the exponent of the error bound for 

{x G M"| P{x) < 0} at the origin is at most r = ■ Moreover, there is a constant cq > 0 such 

that for all x around x = 0 we have the relationships 


Amax(^’(a;)) = max ■(max' 


x^'^ - X,i 


(l+^+ 


2 < i < n 


= max < max ■ 


p 






xf"* - Xi i 

, xf 

+ 

> 

2 < i < n 

> 

(l+Xi- 

L 2 

0 


> 


max I max |co [xf‘^ — Xi_i]_|_, xf| | 2 < i < n|. 


where the second equality is due to xf > 0 while the last inequality follows as the denominator in the third 
equality approaches one as x —> 0. Considering next the finite convex polynomial system 


f < 0 and Co (x"“ — Xi_i) <0 for i = 2,... ,i 


note that its solution set is the same as SpMi = {0}- It is known [S] that for any finite convex polynomial 
system, where all the involved polynomials have degree at most do on K”, a local Holder error bound holds 
with the exponent ( 2 do-iy+i ■ H^nce there are nnmbers c, 5 > 0 such that 


d(x, SpMi) < cmax { max {cq [: 


xf^' - Xi 


i] + , xf}\ 2<i<n} 


(4d-l)'‘+l 


^ cAmax(^(^)) 


(4d-l)"'+l 


< <5. 
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Thus the true exponent of the Holder error bound at the origin in this example lies within the interval 
[ ( 4 d-i)"+i ’ d( 2 d)"-i- ]' other hand, the exponent value calculated in Corollary 14.61 is 

1 1 

^ ~ ^{2n'^ + 4n, d + 3) ~ {d + 3)(3d + 6)2"^+4"-i ’ 
which therefore gives us only a lower estimate of the true Holder exponent in the local error bound. 

We now list two specific structures of polynomial systems, where error bounds can be sharpened: 

• Let (x, y) I— fi{x, y) be a concave polynomial for all Z = 1,..., L, 17 be a compact convex set, and there 
exist xq € M" with max fi{xo,y) < 0. It can be directly verified in this case that x i—>■ max fi{x,y) is a 

1<Z<L 

continuous concave function. Note that a Lipschitz local error bound (i.e., with the best possible exponent 
T = 1) holds for a concave function satisfying the imposed strict feasibility condition; see jSS] Theorem 3]. 

• Let X I—>■ fi{x,y) be a convex polynomial for each y G K™ and I = 1,... ,L, and let 17 be a convex 

polytope in M™. Since in this case the number of extreme points of 17 is finite and is at most to + 1, 
we get the representation 17 = co{y^^\... Noting that a convex function attains its maximum 

over a convex polytope at its extreme points, it follows that max f{x,y) = max fi{x,y^'^'^). Thus the 

l<i<m+l 

1<1<L — — 

corresponding explicit error bounds can be simplified and sharpened in this case by using the results from 
our preceding paper for finite index sets 17. 

5 Higher-Order Stability Analysis for Optimization Problems 

In this section we develop some applications of the error bound results from Section 4 to higher-order 
stability analysis for three important classes in parametric polynomial optimization dealing with infinite 
sets of parameters Y{x) from (11.61) . In contrast to the majority of publications in parametric optimization, 
our main emphasis is on higher-order stability of set-valued solution maps with explicit calculating the 
corresponding exponents. Let us start with problems of generalized semi-infinite programming. 

5.1 Generalized Semi-Infinite Programs with Polynomial Data 

Consider the following polynomial generalized semi-infinite program: 

(GSIP) minimizea;gRn poix), 

subject to fi{x, y) <0 for all y G Y{x), I = 1,..., L, 

where the sets Y{x) C are taken from (11.61) . and where all the functions po, fi, gi, and hj are real 
polynomials with degrees at most d. Along with problem (GSIP) we consider also its perturbation 

(GSIPu) minimizea;gRn poix,u), 

subject to fi{x,y)<0 for all y £ Y{x), I = 1,..., L, 

where the cost function p(x,u) with p(x,0) = po(x) depends on the perturbation parameter u G being 
polynomial of degree at most d in x and locally Lipschitzian in u. Note that the case oip{x, u) = Pq{x) — {u, x) 
and Y{x) = {0} in (GSIPu) corresponds to the so-called tilt perturbation introduced by Poliquin and 
Rockafellar in |33] . where the focus is on single-valuedness and Lipschitz continuity of the optimal solution 
map around the point in question, which is known as tilt stability. In mus] the reader can find recent 
second-order characterizations of tilt stability in the aforementioned setting of nonlinear programming with 
C^-smooth data and related issues. The quantitative higher-order stability for polynomial GSIPs established 
in the next theorem is significantly different: we justify a certain Holder continuity of the set-valued solution 
map at the reference point with calculating the corresponding Holder exponent. 
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Theorem 5.1 (quantitative higher-order stability for polynomial GSIPs). Let the feasible set 

G := {a; G ]R"| //(x, y) < 0 for all y GY {x), Z = (5-1) 

be nonempty and compact. Assume that the mapping Y in (I53D is i.s.c. on K" with nonempty and uniformly 
bounded values and that MMFCQ holds on R” relative to Y{x). Denote by S(u) the set of optimal solutions 
to the perturbed problem {GSIPu). Then for any fixed u there are constants c,e>0 such that 

S{u) C S{u) + c ||m — m||'^ if ||u — u|| < e with r = A?(2n + (m + r + s + 2){n + 1), d + L + 2)~^. (5.2) 


Proof. Consider the optimal value function in {GSIPu) given by 

p^,{u) := ini {p{x,u)\ fi{x,y) < 0 whenever y G Y{x), I = 1,..., L}, uGR^, 
and observe that the solution map S{u) to (GSIPu) is represented as 

S{u) = {x G M" I fi{x^ y) F 0 for all y G Y{x), I = I,... ^ and p{x, u) — p.t{u) < 0}. (5-3) 

Denote fi+i{x, y) := p(x, u) — Pt:{u) and define 

fix) := ma.-K^p{x,u) - p^{u), ^max /i(x,y)|. (5.4) 

Applying the error bound from the second part of Theorem 14.21 to the system in (15.31) and (15.41) with L 
replaced by L -I- 1 > 2, we find cq > 0 such that 

dist(a;, S'(? 2 )) < Co [(/('(a;)] ^ for all xGG with t = t%(2n + {m + r + s + 2){n + l),d + L + 2) 

Since p{x, •) is locally Lipschitzian around u and G is compact, we find constants e > 0 such that 

|p(a:, ui) — p(a;, M 2 )| < f||Mi — U 2 II for all xGG and ui,U 2 GM^^u). (5.5) 


Denote c := (2/3“^£)'^ with /3 := Cg ^ >0 and for any w G S{u) select z G S{u) with ||w —z|| = dist(w, S'(u)). 
To justify the Holder stability (j5.2l) . it remains to verify the estimate 

||w; — z|| < c ||u — ■u|l'^. (5.6) 

We consider the case of ic G S{u)\S{u)-, otherwise, (15.61) is trivial. Then 

p{w,u) -p^,{u) = (j){w) = [4>{w)]^ > Pd{w,S{u))'^ = P\\w - z\\i. 

Since z G S{u), we get p{z,u) = p^:{u) < p{w,u), and hence 

\\w - z\\i < {p{w,u) - p4u)) < f~^{p{w,u) - p{z,u)). (5.7) 

Furthermore, it follows from w G S{u) that p{w,u) < p(z,u), and so (15.5|) gives us the relationships 

p{w,u) - p{z,u) = {p{w,u)-p{z,u)) + {p{z,u)-p{z,u))+ {p{w,u)-p{w,u)) 

< {p{z , u) — p{z, u)) + {p{w, u) — p{w, m)) < 2i\\u — u|| as w,z G G 

implying together with (15.71) that Hw — z\\i < j3~^{p{w, u) — p(z, u)) < 2j3~^l\\u — m||. Therefore 

dist(r(;, S{u)) = Hw — z|| < c \\u — u\\'^, 

which justifies (j5.6|) and thus completes the proof of the theorem. □ 

Remark 5.2 (sharper exponents in higher-order stability for equality descriptions.) Based on 
Remark 01 observe from the proof of Theorem 15.11 that the Holder stability exponent in (OD can be 
sharpened to t = ^{2n + {m + r + s + 2){n + 1), d + L + 1) if the sets Y(x) in {GSIP) are described by the 
polynomial equalities only: Y(x) = {y\ hj {x, y) = 0 , j = 1 ,..., s}. 


17 


5.2 Optimization of Polynomial Matrix Inequalities 

Next we address stability issues in optimization of polynomial matrix inequalities: 

{PMI) minimize po{x) subject to P{x) ^ 0, 

where po: R." —>■ R. is a real polynomial of degree not exceeded d, and where the mapping P: R" —>■ 5'™ 
is such that each (*,j)th element of its image {P{x))^^, 1 < i,j < m, is a real polynomial under the same 
degree assumption. As in Subsection 5.1, consider the perturbed version of (PMI): 

{PMIu) minimize p{x,u) subject to P{x) ^ 0, 

where u S R'^, p{x, 0) = po(a;), the degree of p{-, u) is at most d, and p{x, •) is locally Lipschitzian. 

Theorem 5.3 (higher-order stability of optimal solution sets to polynomial matrix inequalities). 

Assuming that the feasible set G := {a:| P{x) ^ 0} in {PMIy) is nonempty and compact and then denoting 
by S{u) the set of optimal solutions to this problem, for each u we find constants c,e > 0 such that 

S{u) C S{u) + c ||u — 'u||'^ if ||u — ull < e with r = i%{2n + (m + 3)(n + 1), d + 4) 

Proof. Consider the function f{x, y) := Vi which is a polynomial of degrees at most d + 2. 

We clearly have the equivalence 

P{x) P 0 f{x,y) < 0 for all \\y\\^ = 1. 

Defining now the constant sets Y (x) by 

r(x) = D:={yeR-| \\yf = l}, 

observe that the set D is regular. Then applying Theorem 15.11 with L = 1, r = 0, and s = 1, replacing d 
with d + 2, and taking into account Remark 15.21 justify the claimed Holder stability. □ 


5.3 Second-Order Cone Programs with Polynomial Data 


The class of polynomial optimization problems considers below belongs to second-order cone programming 

{SOCP) minimize Po (a;) subject to \\{f[‘\x),..., f^\x))\\ < f^\j^{x), l = l,...,L, 

where po, fj^'^: R" —^ R for j = 1,..., m + 1 and I = 1,... ,L are real polynomials with degrees at most d. 
Besides other areas of applications, problems of this type arise in optimization under uncertainty via the 
robust optimization approach. As a simple illustration, let us consider the following quadratic optimization 
problem under uncertainty in constraint data: 

minimize x'^Ax + a^x subject to x'^Bix+ bjx + Pi <t), l = l,...,L, (5.8) 


where the triples {Bi, bi, Pi) G S'" x R" x M., I = 1,..., L are not known in advance. As commonly accepted, 
suppose that the constraint data belong to the ellipsoidal uncertainty set 


Ui := 


{{Bf,bf\pf^)+Y^u\{B\ 






The robust optimization approach deals with an associated problem of minimizing an 
uncertainty. In our case it reads as follows: 


adverse effect of 


minimize xF Ax + a^ x subject to x^ Bix + bf x + Pi < 0 for all {Bi,bi, Pi) gUi, l = l,...,L. 
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By direct calculations we have for each I = 1,... ,L that 


sup {x^Bix + bf X + Pi} 
{Bi, bi, Pi) S Hi 


= x"’" + {b[°'^)'^x + pI^'^ + sup v}(x^b}^^x + {b\^'^)'^ X + p\^^) 

IIK.-.«f)ll<ij=i 

= x^b}°^x + {b\°^)'^x + p}^'^ + II }x^b}^^x + + pI^\ . .., x^b}^^x + {b\^'^)'^x + ||- 

This therefore reduces the uncertainly model to the above framework of (SOCP). 

Now we consider the following perturbations of {SOCP) given for all m S R'^ by 

{SOCPu) minimize p(x,m)) subject to || (/f^(x),...,/£^(a:)) || </®+i(a;), 1 = 




,(o) 


where p{x, 0) = po( 2 ;), p{-,u) is a polynomial of degree at most d, and p{x, ■) is locally Lipschitzian on R'J. 

Theorem 5.4 (higher-order stability of optimal solution sets to polynomial second-order cone 
programs). Assuming that the feasible set G = {x\ \\{fi\x),...,fm{x))\\ < f^+ii^)} {SOCPu) is 
nonempty and compact, for any u G denote by S{u) the set of optimal solutions to this problem. Then 
given u gW^, there exist constants c, e > 0 such that 

S'(m) C S')?!) + c ||u — mII"^ if ||u —ull^e with t = I%{2n + {mL + L + 2){n + l),d + L + 1) (5.9) 


Proof. For each I = 1,..., L we define the function 


fi(.x, y) ■■= yjfj{x) - f^\i{x), 
i=i 

which is a polynomial with real coefficients on M"+"* of degree at most d + 1. Note that 

m 

\\{fi\x),...,f}}}{x)) \ < f^\i{x) ^ sup|y]]yj'Vi'^(a;) - fL+iix) || (yf \ , y^^) || = l} < 0 

3 = 1 

whenever I = 1,... ,L. Consider the sets Y {x) C R™^ given by 

Y{x) =n:={y= (yW,... ,y(^))| = {y[^\ ..., y«), = 1, I = 1,..., L}, x G RC 


Observe that these sets are nonempty and uniformly bounded around any x G R" and that the set O is 
regular. Applying now Theorem 15.11 and Bemark 15.21 with r = 0 and s = L and with replacing m by mL, 
we arrive at the claimed Holder stability in (15.9p . □ 


6 Convergence Rates in Feasibility and Asymptotic Analysis 

The concluding section of the paper presents applications of the obtained error bounds to deriving explicit 
convergence rate estimates in two important frameworks involving polynomial functions. The first issue 
concerns the cyclic projection algorithm to solve the well-known convex feasibility problem, while the second 
one deals with asymptotic analysis of subgradient dynamical systems. 

6.1 Cyclic Projection Method for Convex Feasibility Problems 

A common problem arising in diverse areas of mathematics and applications is to find a point in the in¬ 
tersection of closed convex sets Ci, I = 1,..., L. This problem is often referred to as the convex feasibility 
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problem. One of the most powerful methods for solving the convex feasibility problem is the so-called cyclic 
projection algorithm, which is formulated as follows. Given finitely many closed convex sets Ci,.. .Cl in M" 
with nz=i C'z ^ 0 , pick xq S R" and denote tt; := nci for 1 = 1 , • ‘ ‘ jL, where nci stands for the (unique) 
Euclidean projection to the set Ci. The sequence of cyclic projections {xk)k&i is defined by 

xi := TTiXo, ...,XL ■■= ttlXl-1, xl+i ■= ttiXl, - ( 6 . 1 ) 


In the case of L = 2 the cyclic projection method reduces to the well-known von Neumann alternating 
projection method] see, e.g., 011] and the references therein, where the reader can find various practical 
applications of the cyclic projection method and related algorithms. 

Convergence properties of the cyclic projection algorithm have been examined by many researchers. In 
particular, the seminal paper by Bregman 0 shows that the sequence {xk)k&i generated by cyclic projections 
always converges to a point in C. Furthermore, it is proved in m that the cyclic projection algorithm 
converges linearly for regular intersections of convex sets. On the other hand, for convex sets with irregular 
intersections (e.g., when the intersection is a singleton), the cyclic projection algorithm may not exhibit 
linear convergence even in simple two-dimensional cases as observed by 0 Example 5.3]. This raises the 
following basic question: Can we estimate the convergence rate of the cyclic projection algorithm for convex 
sets with possibly irregular intersections of certain particular structures? Quite recently 0 this issue has 
been investigated in the setting where the convex sets Ci are the so-called basic semialgebraic given by 

C'i:={a:GR"|g](a:)<0, i = l,...,n}, 1 = 1,...,L, (6.2) 

with convex polynomials gl of degrees at most d. It is shown in 0 that for d > 1 the sequence of cyclic 
projections {xk)ken from (ED) converges at the rate of 


^ min|(2d 1)” I, 2(j^^” 

where the symbol [a] stands for the integer part of the number a. 

In what follows we establish an explicit convergence rate estimate for the more general case of convex 
sets described by polynomial matrix inequalities. It is said that a set C is polynomial matrix inequality 
representable if there are numbers m,d G N and a mapping P: R" —>■ S'™ such that each (i,j)th. element of 
its image (P(a:))^^., I < z, j < m is a real polynomial with degree at most d and we have the representation 
C = {a: G R"! PIx) ^0}. It is obvious that any basic semialgebraic set (16.21) belongs to this class. The 
following example shows that the converse is not true. 

Example 6.1 (polynomial matrix inequality representable sets may not be basic semialgebraic). 

Given an odd number d, consider the set 


C := 


|(xi,a: 2 ) G R^ 


-xf 1 

1 -xi 



= {ixi,X2)\ X 1 X 2 > 1 , 


Xl >Q,X2> 0 }, 


which is clearly convex and polynomial matrix inequality representable. Let us show that C cannot be 
written in the basic semialgebraic form H6.2I) . Assuming the contrary, take the function f(xi,X 2 ) ■= xi 
for which f{x) = 0. It follows from the polynomial extension of the celebrated (quadratic) 

Frank-Wolfe theorem that any convex polynomial / bounded from below over a basic semialgebraic convex 
set attains its minimum; see 0. This gives us a: G C such that f{x) = f{x) = 0 for the function / 

and the set C defined above. It surely contradicts the structures of / and C, and thus verifies the claim. 

To derive an explicit convergence rate estimate for the cyclic projection algorithm applied to polynomial 
matrix inequality representable convex sets, we use the error bound result from Theorem l4.2l and the following 
abstract result on cyclic projections established in 0 under a certain Holder condition on convex sets. 
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Lemma 6.2 (conditional cyclic convergence rate). Let Ci C K", I = he closed convex sets 

with C := C'i 7 ^ 0, and let {xk)k^n be the sequence of cyclic projections (16.11) with a starting point Xq. 
Impose the following Holder condition of rank r G (0,1]; for any compact set K C M" there is a positive 
number cq providing the estimate 

di{x,C) < co( dist(a:, Cj) I whenever x G K. (6-3) 

1=1 ' 

Then there exists Xoo G C such that the sequence (xk)kGn converges to Xoo, and for some M > 0 we have 

1 

\\xk - Xoc\\ < M k 2t-1-2, 


Now we are now ready to obtain an explicit convergence rate estimate of the cyclic projection algorithm 
for polynomial matrix inequality representable convex sets. 

Theorem 6.3 (explicit convergence rate of the cyclic projection algorithm). Let Ci be polynomial 
matrix inequality representable convex sets given in the form 

a := {x G M”| A®(a;) ^ 0}, Z = 1,...,L, 

where every ^ ^ matrix mapping such that each 1 < i,j < m, is a real polynomial 

with degree at most d > 1. Given xg G K”, consider the sequence (xfc)fcgN of cyclic projections (16.11) . Then 
Xk converges to some x^o G C and there is M > 0 ensuring the estimate 

\\xk - a;oo|| < Af k gN, 

where p := [2 3£{2n + (to + 3)(n + l),d + L) — 2\~^ with the constant M defined in (12.41) . 


Proof. Let us first show that for any compact set K C M" there is cq > 0 such that condition (16.311 holds 
with some explicitly calculated r G (0,1). To proceed, fix a compact set K and define 


n := {y G M"*! IIpII^ = 1} and fi{x,y) := ^ y^ (v4(')(x)) yj, l = l,...,L. 

i,l=l 

Then we readily have the polynomial set representations 

C; = {x G fi{x,y) < 0 for all y G Ll} and C = |x G M” 

1<1<L 

Applying now the error bound result of Theorem 14.21 to the set C from (16.41) with Y{x) = ft, r = 0, and 
s = 1 while taking into account Remark 14.31 and the regularity of the set LI above, we find cq such that 


max /i(x,y) < o|. (6.4) 


dist(x, C) < Cq 


(r 

^ 1<1 


max 

<L 




< Co 


L 

E 


[max/i(x,y)] 

ygn 


for all X G K, 


where r = .^(2n + (to + 3)(n + l),d + L)~^. For each I = 1,...,L define the function hi{x) := 
[maxy^Q fi{x,y)] + , x G R", and observe that it is locally Lipschitzian by Proposition 13.11 Due to the 
compactness of K we have Mq := max^^K d{x, Ci) < oo. Since the set K := K + MqE is also compact, 
we get c; > 0 ensuring the global Lipschitz condition on this set: 

|/i/(xi) -/i;(x2)| < Ci||xi - X2II for all X1,X2 G AT, / = 1 ,...,L. 
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Now pick an arbitrary vector x £ K C K and select xi £ Ci such that dist(a;,C';) = ||x — Xi||. Denoting 
c := maxjci,..., cl} > 0 and observing that Xi £ K and hi{xi) = 0 give us the relationships 


L L 

^dist(a;,Ci) = ^ \\x-xi\\ 
1=1 1=1 

which ensure therefore the estimates 

dist(a;,C') < cq 


> 

1=1 ^ 1=1 


L \ T 




h{x) j < Co ( c dist(x, Cl) 


1^1 




Thus we arrive at (|6.3I) and complete the proof of the theorem by using Lemma 16.21 


□ 


6.2 Subgradient Dynamical Systems 

The final piece of this paper is devoted to applications of the obtained error bounds to conducting a quan¬ 
titative asymptotic analysis of the subgradient dynamical systems given by 

. . Q £ x{t) + d(t){x{t)) a.e. on (0,T), 

d(j)[x{t)) 7 ^ 0 for all t £ [0,T), 

where </> is a semialgebraic function on R" in the form (j){x) := sup{/(a;,?/)| y £ D}. We suppose in what 
follows that D is a regular, compact, and semialgebraic set defined by 

D := {y e R™| g^{y) <0 for i = 1,... ,r for hj{y) =0, j = 1,... ,s} (6.5) 

and that all the functions /, gi, and hj are real polynomials with degrees at most d. By a solution of 
(DS) we understand any absolutely continuous curve x: [0, T) —>• R" satisfying the relationships therein. A 
trajectory x{t) is called maximal if there is no possible extensions of its domain compatible with {DS). Our 
asymptotic analysis of trajectories for {DS) as t —>■ oo is based on Theorem 13.81 extending the Lojasiewicz 
inequality to maximum functions and on the following result taken from [7]. 

Lemma 6.4 (convergence to critical points). Let x{t) be a bounded maximal trajectory of{DS), where 
the function (j) is assumed to be bounded from below on K". Then x{t) is defined on [0, oo) and converges to 
a critical point x £ R." of (j), i.e., such that 0 £ d(j){x). Take any 9 £ [0,1) for which the quantity 
is bounded for any x around x under the convention that ^ = 1. Then there are numbers c, to >0 such that 

||a;(t) — 5|| < c (t-I- for all t > to- 

We finally arrive at the constructive rate evaluation of asymptotic convergence for maximal trajectories. 

Theorem 6.5 (explicit rate estimate of asymptotic convergence). Let (j) be a semialgebraic function 
on R” given by (f>{x) := max{/(x, y)| y £ D}, where LI is a semialgebraic set in form (16.51) under the 
assumptions formulated above. Suppose in addition that cf is bounded from below on R". Then any bounded 
maximal trajectory x{t) for {DS) converges as t ^ oo to a critical point x £ R" of f). Furthermore, there 
are numbers c, to > 0, which ensure the estimates 

||x(t) — x|| < c (t -I- and \(l){x{t)) — (j){x)\ < ct~^^ for all t > to, (6.6) 

where 9 = 1— i%{2n {m + r s){n -I- 1), d -I- 2)“^ with the constant Si defined in (12.411 . 
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Proof. Let a; be a critical point of (f) the existence of which follows from Lemma [6.41 Denote r := ^{2n + 
{m + r + s){n + 1), d + 2)~^ and 9 := 1 — t. Then Theorem l3. 81 tells us that there are numbers cq, £ > 0 such 
that m0(a::) > c|^(a;) — 4>{x)f for all x G ]Be(a;), and therefore 


{(fix) - (j){x )\'^ 
m^{x) 


< Co 


whenever x G Be(a:). 


(6.7) 


Thus the first estimate in () 6 . 6|1 follows directly from Lemma 16.41 Define now 7 : [0, 00 ) —^ R by 


'y{t) := (j>{x{t)) - (j){x), tG[0,oo), 


and observe that both functions (f and 7 are Lipschitz continuous and so a.e. differentiable on (0,oo). This 
yields d4>{x(t)) = {V(^(<)} for a.e. t G (0,oo); see, e.g., [27l Corollary 1.82]. Then we see that 

i{t) = {'^<i){x{t)),x(t)) = -||V<()(a:(t))||^ = -m^{x{t)f < -cl\(j){x{t)) - (j){x)\^^ = -cg7(t)^® < 0, (6.8) 

where the second equality follows from 0 G x{t) + d(j){x{t)) while the first inequality follows from (16.711 . This 
shows, in particular, that the function 'yft) is nonincreasing. Since x(t) —>■ x and 7(f) —>■ 0 as f —>■ 00, we get 
that either j{t) = 0, or j{t) > 0 for all large t > 0. Integrating (|6.8I) allows us to find c, to >0 such that the 
second estimate in (16.61) holds for all t > to, and thus we complete the proof of the theorem. □ 
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